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1. Introduction 

An uncharted area of physical and observational cosmology is the physics of ultra-large scales. 
By this we mean length scales which are near or beyond the cosmic horizon. Most large-scale 
surveys have been limited to sampling length scales of order hundreds of megaparsecs, whilst 
deeper surveys (i.e. higher redshift surveys) are not yet good enough to accurately measure modes 
which are of order the horizon. Yet, the cosmological information hidden in the extremely large 
scales is of great importance for our comprehension of the Universe. For a start, the growth of 
cosmological perturbations is well within the linear regime on those ultra-large scales. This allows 
for better checks of our theoretical model, since we can there safely disregard the non-linear growth 
of structures—which always requires some degree of ad hoc modelling. Furthermore, on very large 
scales the most relevant physics is simply given by the gravitational interaction, thus freeing us from 
the necessity of taking astrophysical processes into account. 

Theoretical knowledge of the Universe on the largest cosmic scales can be very precise, but 
ultra-large scales have been most difficult to access. This is mainly due to two observational issues. 
Firstly, if we want to probe very large scales—i.e. very large wavelengths or very small Fourier 
momenta—a wide field of view is nof sufficienf. For a wide-held, shallow survey fhe fransverse 
componenfs, kj^, of fhe physical wavevecfor k may be arbifrarily small, buf fhe smallesf measured 
wavenumber, k = |k| = y/k^ -|- is limifed by fhe facf fhaf fhe survey probes only a fhin shell— 
fhus providing us wifh a large parallel componenf, k\\. In ofher words, we need fo observe a very 
large ‘cube’ of Universe. This is a challenging fask for convenfional galaxy redshiff surveys, for 
if is exfremely hard fo achieve fhe required sensifivify af high redshiff over a large area of fhe sky. 
Closely relafed fo fhis is fhe second observafional problem, namely cosmic variance. Even fhough 
fhe full SKA will be able fo defecf HI galaxies as far away as z ~ 2 over roughly fhree quarfers of 
fhe celesfial dome—fhus acfually observing a huge cosmic cube—fhe relafively small number of 
sources wifh ulfra-large separafion can heavily limif fhe consfraining pofenfial of fhe survey. 

The SKA is in a unique posifion fo change fhis sfafe of affairs and push forward fhe sfudy 
of ulfra-large scale modes. To fackle successfully fhe fwo observafional issues described above, 
fhe SKA will employ fwo newly developed fechniques, one being HI infensify mapping (IM) and 
fhe ofher fhe so-called mulfi-fracer approach. Since we can exploif IM bofh during and affer fhe 
epoch of reionisafion (EoR), fhe SKA will be able fo probe excepfionally large redshiffs. When 
mapping large-scale sfrucfure affer fhe EoR, if is possible fo efficienfly collecf informafion on large 
wavelengfhs wifhouf having fo worry abouf resolving small-scale feafures or individual galaxies. 
When mapping fhe EoR, one can effectively sample a large number of horizon-size volumes af fhe 
time. 

In addifion, fhe mulfi-fracer fechnique will allow us fo beaf down cosmic variance. This ap¬ 
proach is based on galaxies being biased buf nof stochastic fracers of fhe underlying densify field. 
Thus, by observing differenf galaxy populafions—namely, differenf fracers—we can consfrucf sev¬ 
eral realisations of fhe clusfered halo disfribufion. Alfhough measuremenfs of global cosmological 
properties like fhe fofal matter densify Q.m or fhe Hubble consfanf Hq will sfill be affecfed by cosmic 
variance, we will be able fo consfrain halo properfies to a much higher degree of accuracy. 

In fhis Chapfer, we oufline fhe science fo be done by using fhe SKA fo map fhe matter disfribu¬ 
fion on fhe largesf cosmic scales. We firsf make fhe case for why ulfra-large scales are inferesfing. 
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focussing on relativistic effects, modified gravity (MG) and scale-dependent biasing arising from 
primordial non-Gaussianity (PNG). Then, we address how different methods will constrain these 
various physical phenomena, showing that indeed the SKA will be transformational for constrain¬ 
ing PNG, and unique for probing gravity on super-horizon scales. Lastly, we discuss the issues that 
need to be faced when trying to probe the largest scales. We refer the reader to Chapters Santos 
et al. (2014), Abdalla et al. (2014) and Jarvis et al. (2014) for details on the phases of the SKA, its 
observation modes and the different methods we are discussing here. 

2. The Importance of Probing Ultra-Large Cosmic Scales 

The study of the biggest volumes ever of cosmic large-scale structure will allow a major ad¬ 
vance in tackling two of the most fundamental questions in contemporary cosmology: What is the 
physical mechanism that generated the inflationary expansion in the early Universe? Does general 
relativity (GR) hold on the largest scales? 

Most models of inflation do predict some degree of non-Gaussianity in the distribution of 
primordial density fluctuations, but the simplest models predict negligible PNG. A detection of 
PNG will be vital for ruling out classes of models and advancing our understanding of inflation. 
Currently, the most stringent constraints on PNG come from the Planck satellite, but future cosmic 
microwave background (CMB) experiments are unlikely to improve these results significantly. The 
new frontier of constraining PNG is large-volume surveys of the matter distribution. 

Similarly, tests of GR on cosmological scales are based on observations of the large-scale 
structure. Current constraints are weak, but with its huge surveyed volumes and multiple probes, 
the SKA will take the lead in the next generation of tests. Additionally, we can tighten the current 
constraints on dark energy (DE) and MG models by including much larger scales, thus increas¬ 
ing the statistical power of the observations and improving the sensitivity to any scale-dependent 
deviation from GR. 

In this Section we review such ultra-large scale phenomena, highlighting their importance for 
our understanding of the Universe. 

2.1 General Relativistic Effects 

Probing ultra-large scales involves a theoretical challenge that has been recognised only re¬ 
cently. GR effects arise from observing on the past light-cone, which distorts the number counts 
and brightness temperature fluctuations on very large scales. So far, most analyses have been 
performed using a Newtonian approximation, with the flat-sky redshift-space distortions (RSDs) 
grafted on as a special relativistic effect. This is the ‘standard’ relativistic correction to the New¬ 
tonian approximation, which is significant on sub-horizon scales. Sometimes, the effect of lensing 
convergence is also included. This is adequate for past and present surveys, which analyse galaxy 
clustering on scales well below the horizon. Sometimes this Newtonian-like approach also includes 
the contribution of weak lensing magnification to the matter over-density. The lensing effects can 
be significant on small scales. But future wide and deep surveys will need to employ a more pre¬ 
cise modelling, accounting not only for RSDs and lensing, but for all geometric and relativistic 
corrections. The full relativistic analysis includes terms that are suppressed on sub-horizon scales 
such as velocity (or Doppler) terms, Sachs-Wolfe (SW) and integrated SW (ISW) type terms, and 
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time-delay contributions (e.g. Yoo 2010; Jeong et al. 2012; Bonvin & Durrer 2011; Challinor & 
Lewis 2011; Bruni et al. 2012). 

One can work with the redshift-space correlation functions and include all GR and wide-angle 
contributions (Bertacca et al. 2012; Raccanelli et al. 2014a, 2013; Bonvin et al. 2014). By analysing 
this, we can potentially extract more information about the structure of galaxy clustering. For 
instance, the odd Legendre multipoles of the correlation function vanish in the Newtonian flat-sky 
approximation, whereas they are in general non-zero for wide-angle separations and in the full GR 
analysis (Raccanelli et al. 2014a; Bonvin et al. 2014). Alternatively, the use of the angular power 
spectrum also avoids the flat-sky assumption, which is important for full-sky surveys with many 
wide-angle correlations (Bonvin & Durrer 2011; Challinor & Lewis 2011; Di Dio et al. 2013). 

To give a flavour of how GR corrections alter the Newtonian prediction, we show in Fig. |] 
the ratio between a galaxy clustering angular power spectrum in which some GR effect has been 
switched on and the simple Newtonian result. We factorise the various corrections as follows: 
RSDs (top, left panel); velocity (top, right panel), where we include terms proportional to the 
velocity along the line of sight; weak lensing convergence (bottom, left panel); and potentials 
(bottom, right panel), which account for the effect of gravitational potentials at the source, SW, ISW 
and time-delay terms (see e.g. Challinor & Lewis 2011). The two sets of curves refer to different 
surveys, all of them assumed to observe a certain population of galaxies with constant bias 1.5, 
distributed according to a Gaussian window centred at redshift Zm and with width o^- Specifically, 
we compare a shallow survey with Zm = 0.5 (red curves) to a deep survey with Zm = 2.0 (blue 
curves). Moreover, we consider both the cases of a narrow (solid lines) and a broad (dashed lines) 
window, viz. = 0.01 and 0.1 for the shallow survey and Ow = 0.05 and 0.5 for the deep survey. 
From Fig. [T], it is apparent that, although RSDs and lensing are the dominant, the other terms can 
also become important, particularly for as deep a survey as the SKA. 

The GR corrections to the power spectrum should be observable for certain SKA surveys, 
using the multi-tracer technique—the task of forecasting this capacity is under way. Even if these 
corrections are not directly detectable, it is essential to include them in theoretical analysis of the 
power spectrum, for two related reasons: 

1. to avoid biasing parameter estimations through a theoretical systematic of incorrect mod¬ 
elling; 

2. to correctly extract maximal information from the largest scales. 

2.2 Primordial Non-Gaussianity 

The largest cosmic scales are a crucial source of information about the physical processes at 
play during the early stages of the Universe’s evolution. The standard model of inflation and its 
generalisations predict seed primordial density fluctuations with some level of non-Gaussianity in 
the probability distribution (e.g. Bartolo et al. 2004). We can parameterise the non-Gaussianity 
in Bardeen’s gauge invariant potential <I> as the sum of a linear Gaussian term and a quadratic 
correction (Verde et al. 2000; Komatsu & Spergel 2001), i.e. 

<I> = 0+/nl(0"-(0')). (2.1) 
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RSD Velocity 





Figure 1: Ratios of angular power spectra including some kind of GR corrections to the standard Newtonian 
approximation, Red (blue) curves refer to a shallow (deep) survey, while solid (dashed) lines to 

natTow (broad) source redshift distributions. (See text for details.) 


If the distribution of primordial density perturbations is not Gaussian, it cannot be fully described 
by a power spectrum; we rather need higher-order moments such as the bispectrum. In particular, 
different models of inflation give rise to different bispectrum shapes, thus making the study of PNG 
valuable for obtaining a deeper knowledge of the physics of inflation. 

The standard single-field inflationary scenario generates negligibly small deviations from Gaus- 
sianity. These deviations are said to be of the local shape, and the related bispectrum of Bardeen’s 
potential is maximised for squeezed configurations, where one of the three wavenumbers has a 
much smaller magnitude than the other two. In this case, the PNG parameter, /nl, is expected 
to be of the same order as the slow-roll parameters, namely very close to zero (Falk et al. 1993). 
However, this does not mean that /nl ~ 0. Due to the inherent non-linearity of GR, it is not pos¬ 
sible for Gaussianity in the primordial curvature perturbation to be reflected exactly in the density 
perturbation. A non-linear GR correction to the initial conditions leads to an effective = —5/3 

in large-scale structure (Verde & Matarrese 2009). 

Local-shape PNG can also be generated when an additional light scalar field other than the 
inflaton contributes to the observed curvature perturbations (Bartolo et al. 2004). This happens for 
instance in curvaton models (Sasaki et al. 2006; Assadullahi et al. 2007) or in multi-field mod¬ 
els (Bartolo et al. 2002; Bernardeau & Uzan 2002; Huang 2009). Other than local-type PNG, 
there are inflationary models in which the kinetic term of the inflaton Lagrangian is non-standard. 
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containing higher-order derivatives of the field itself. One significant example of this is the ‘DBF 
model (Alishahiha et al. 2004; Arkani-Hamed et al. 2004; Seery & Lidsey 2005), where the primor¬ 
dial bispectrum is maximised for configurations where the three wavevectors have approximately 
the same amplitude—the so-called equilateral-type PNG (Creminelli et al. 2007). Otherwise, for 
deviations from Gaussianity evaluated in the regular Bunch-Davies vacuum state, the primordial 
potential bispectrum is of local or equilateral shape, depending on whether or not higher-order 
derivatives play a significant role in the evolution of the inflaton field. If the Bunch-Davies vacuum 
hypothesis is dropped, the resulting bispectrum is maximal for squashed configurations (Chen et al. 
2007; Holman & Tolley 2008). Lastly, another shape of the bispectrum can be constructed that is 
nearly orthogonal to both the local and equilateral forms (Senatore et al. 2010). 

The best probe of PNG up to now has relied on measuring the CMB temperature anisotropy 
bispectrum (see Ade et al. 2014, for the latest results). However, it has been demonstrated that PNG 
also induces an additional, peculiar scale and redshift dependence in a biased tracer of the under¬ 
lying matter distribution (Dalai et al. 2008; Matarrese & Verde 2008; Schmidt & Kamionkowski 
2010; Desjacques et al. 2011). The modification Abx{z,k) to the Gaussian large-scale bias bx of a 
biased tracer X induced by local-type PNG is 

Abx{z,k) = bx{z) + 3[bx{z) - (2.2) 

where bx{z) is assumed scale-independent, dc is the critical value of the matter over-density at 
collapse, the transfer function T (k) —)■ 1 on large scales, and D+ (z) is the linear growth factor of 
density perturbations normalised to unity today. ^ 

It is clear from this equation that the effect of PNG on the power spectrum grows on large 
scales, as k^^. These are the same scales on which GR corrections are becoming significant. 
Therefore, one needs to incorporate the GR corrections in theoretical analysis in order to make 
accurate (and non-biased) predictions and estimates of PNG (Camera et al. 2014b). Similarly, 
an incorrect treatment of PNG on scales where its effects are not negligible may threaten future 
cosmological experiment accuracy (e.g. Camera et al. 2014a). 

2.3 Modified Gravity 

Einstein’s theory of GR has been tested to exquisite precision on ‘small’ scales, namely in the 
laboratory, the Solar System and with the help of pulsars. On the contrary, tests on scales approach¬ 
ing the cosmological horizon are still rather weak. However, there is great interest in testing GR 
on very large scales in the context of DE, as the cosmological constant suffers from a plethora of 
theoretical problems. Since the accelerated expansion of the Universe is a late-time phenomenon, 
it is natural to look for hints concerning the physical nature of the underlying mechanism on very 
large scales. 

The landscape of MG theories is rather heterogeneous, and there is little reason to prefer any 
one model over the others. Attempts to parameterise this landscape in a suitably generic way have 
been only partially successful (Hu & Sawicki 2007; kinder & Cahn 2007; Baker et al. 2013; Battye 

'Note that we adopt the large-scale structure (LSS) convention, which implies ~ (see e.g. Dalai et al. 

2008; Afshordi & Tolley 2008, for more details). 
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& Pearson 2012) although methods based on effeetive field theory approaehes provide a possi¬ 
ble avenue for the large elass of sealar field models with seeond-order equations of motion, like 
Horndeski-type theories (Piazza & Vernizzi 2013; Bellini & Sawieki 2014). Phenomenologieal 
parameterisations ean in prineiple deseribe the full range of deviations from GR relevant for eos- 
mology and thus fully exploit the information eontained in the data by effeetively modelling the 
metrie perturbations (see Kunz 2012, and referenees therein). However, phenomenologieal models 
are mostly useful to eapture the evolution of linear perturbations, whieh effeetively limits them to 
large seales. 

Many MG theories modify the expansion history of the Universe—often by eonstruetion, as 
they are intended as alternatives to DE. However, it has been shown that many DE models are also 
able to produee arbitrary w{z), and MG models ean be tuned to mimie the eosmologieal eonstant 
value of H'(z) = — 1. Therefore, the expansion history is not a smoking gun for MG. On the other 
hand, MG theories in general also alter the growth history of the Universe, i.e. the evolution of 
the metrie perturbations and thus the growth of the matter density eontrast too. Many MG models 
modify the effeetive strength of gravity as a funetion of scale and/or redshift, break the equivalence 
principle and so forth, thus yielding a number of possible observational signatures. On large (lin¬ 
ear) scales, the linear growth rate, f{z), is the most sensitive probe of the growth history, and can 
be well-measured using RSDs from galaxy redshift surveys and IM (see Hall et al. 2013, for IM). 
Another important quantity in this context is the scalar anisotropic stress (or gravitational slip), 
which is generally non-zero in MG theories, especially on very large scales (Saltas & Kunz 2011). 
The anisotropic stress is an attractive test for deviations from the standard model because it is mea¬ 
surable without any assumptions on the initial power spectrum or the bias by combining peculiar 
velocities (e.g. RSDs from IM) with weak lensing measurements (Motta et al. 2013). Its presence 
appears to be also linked directly to a modification of the propagation of gravitational waves, which 
can be used as a way to define whaf MG means (Salfas ef al. 2014; Camera & Nishizawa 2013). 

Going fo very large scales is imporfanf for many reasons. Generally speaking, fhe small scales 
where perfurbafions are highly non-linear lend lo be less useful for cosmological lesls of GR, as MG 
Iheories generically need lo be screened on small scales lo avoid violaling Ihe slrong Solar System 
conslrainls (allhough Ihe mildly non-linear region is inleresling). In addition, accurate predictions 
on non-linear scales require dedicated A-body simulations for each Iheory, and as mentioned above, 
il is unclear how lo use Ihe unifying phenomenological framework on non-linear scales. Baryonic 
effecls on small scales are also badly understood and add a systematic error which effectively 
renders Ihese scales unusable for precision cosmology for Ihe time being. Then, as Ihe effecls lhal 
we are looking for are small, we need as many modes as we can gel, which also pushes us to large 
volumes. Easl bul nol leasl, Ihe scale dependence of deviations from GR is a crucial observable to 
distinguish belween differenl models, and firsl we need a wide range of scales to be able to observe 
a scale dependence, and secondly Ihe horizon scale is a nalural place where to look for Ibis effect 

Evenlually, anolher way to probe MG al early times provides a link to PNG. According to 
Barlolo el al. (2014), modifications of GR during inflation can create a non-zero quasi-local bis- 
peclrum wilh a non-negligible amplilude lhal can be probed wilh Ihe melhods described in Ibis 
Chapter. Eurlhermore, Ihe GR effecls lhal become relevanl on very large scales are aclually benefi¬ 
cial for testing gravity on Ihose scales, as Ihey conlain additional information aboul Ihe evolution of 
perlurbalions. How to exlracl Ibis additional information besl is currenlly an active research topic. 
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but the SKA is uniquely suited to take advantage of this new source of constraints. 


3. Accessing Ultra-Large Cosmic Scales 


In this Section, we review the two envisaged techniques that the SKA will exploit to tackle 
the difficult problem of accessing the largest cosmic scales, viz. IM (both from and after the EoR) 


and galaxy multi-tracing. In Sec. 3.4, we also discuss the most important systematics that we 
shall have to deal with for a fruitful exploitation of such methods. As a figure of merif, we adopf 
fhe achievable accuracy C7(/nl) on a forecasl measuremenf of fhe PNG parameter. However, fhe 
fechniques presented here will in general enable us fo invesfigafe all fhe peculiar phenomena which 
occur on exfremely large cosmic scales. For example, dedicafed analyses aiming af sfudying ofher 
ulfra-large scale effecls such as GR corrections are being performed af fhe fime of fhe edifing of 
fhis Chapfer. 


3.1 Intensity Mapping After the Epoch of Reionisation 

IM is an alternative approach for probing the density field of the large-scale structure (Battye 
et al. 2004; Wyithe & Loeb 2008; Chang et al. 2008; Peterson et al. 2009). It involves mapping out 
the combined emission of the 21 cm, or HI, line from unresolved galaxies. In doing so, the large- 
scale structure is detected in three dimensions. If one foregoes identifying individual galaxies, one 
can greatly speed up the observation and detection of the large-scale structure. IM experiments are 
sensitive to structures at a redshift range that is observationally difficult to span for ground-based 
optical surveys (Seo et al. 2010). SKA IM surveys will cover a large fraction of the sky (around 
25,000 sq. deg.) over an extremely broad redshift range, making it possible to access the largest 
cosmic scales in the late Universe. Here, we outline the basics of late-time IM, but refer to Santos 
et al. (2014) for a more detailed description. 

In order to model the 21 cm power spectrum, it is crucial to quantify the bias function between 
the matter and HI power spectra, which allows us to write 

P^\k,z)=bUk,z)PHk,z), (3.1) 

where P^{k,z) is the underlying matter power spectrum. Several techniques can be adopted in or¬ 
der to extract mock HI power spectra from realistic hydrodynamic or A-body simulations of the 
cosmic large-scale structure. In one of the simplest approaches (Bagla et al. 2010; Guha Sarkar 
et al. 2012; Bharadwaj et al. 2001), the HI is assigned to particles that belong to dark matter haloes 
as identified in the simulated cosmological volume and the HI content of each halo of mass M, 
needs to be characterised. Other, more recent approaches are instead based on a particle- 
by-particle method where the HI is assigned to gas particles according to more refined physical 
prescriptions which take into account self-shielding effects and the conversion to molecular hy¬ 
drogen (e.g. Dave et al. 2013; Popping et al. 2009; Rahmati et al. 2013; Duffy et al. 2012; Mao 
et al. 2012). A comprehensive comparison of several different methods in terms of 21 cm power 
spectrum is redshift space, performed in the post-reionisation era at z ~ 2—4, has been recently pre¬ 
sented by Villaescusa-Navarro et al. (2014), who used high resolution cosmological hydrodynamic 
simulations. They show that the bui{k,z) approaches a constant value to a good approximation at 
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scales k <lh Mpc^', that galactic feedback in the form of winds is not affecting the signal on the 
largest scales and that the signal is dominated at least at z < 6 by HI residing in haloes. Although 
a full characterisation of the bias function in the (k,z)-plane is difficult to make, it is reassuring 
that on extremely large scales the overall HI bias is flat, thus allowing us to perform cosmological 
studies of MG and PNG. It is also worth mentioning that the other key quantity is the total HI 
fraction, flni (see Padmanabhan et al. 2014), which is in fact measured by observations of quasar 
absorption lines of the Lyman-a forest and damped Lyman-a systems at z > 2 (Tescari et al. 2009; 
Zafar et al. 2013). By means of this data we can then constrain the HI bias. 


3.1.1 Primordial Non-Gaussianity Probed by Scale-Dependent Bias 

Camera et al. (2013) explored this method for constraining the PNG of primordial fluctuations. 
IM experiments seem ideally suited for this goal (see Joudaki et al. 2011; Hazra & Sarkar 2012; 
D’Aloisio et al. 2013; Lidz et al. 2013; Mao et al. 2013, for proposals for doing so in the EoR). 
As argued before, the tightest constraints on /nl will be obtained for large and deep surveys. 
Therefore, the volume of the survey determines the ability to probe below /nl of ^(1). For such 
a method to be successful, we need a deep survey with a large bandwidth accessing frequencies 
of 400 MHz and below. Assuming that line-of-sight scattering and self-absorption phenomena can 
be neglected, the HI line radiation can be related to the differential number counts of halo objects 
(e.g. Challinor & Lewis 2011), from which we can estimate bm- Crucially, given our fundamental 
ignorance about the redshift evolution of the bias, we need to span a wide range of redshifts to 
capture bm{k,z) / 1. This is because, as it is clear from Eq. (2.2), no PNG effect is to be detected 
if the tracer we are looking at is unbiased with respect to the underlying matter density distribution. 

Fig. § shows o{f^) contours in the plane of the surveyed area and total observation time 
obtained with an HI IM experiment with bandwidth between 250 and 1000 MHz, corresponding 
to 0.5 < z < 4.5, subdivided into 75 frequency bins of width Av = 10 MHz (Camera et al. 2013). 
Abscissas roughly cover from a 15 x 15deg^ survey to half sky. The three top panels stand for 
the dish survey case, where the y-axis actually shows total observation time, tjOT. multiplied by 
the number of dishes, Nd- Three maximum angular modes are presented, = 25, 60, and 300, 
corresponding to dish diameters of 5, 15 and 80 m at z ~ 3. For higher angular resolution, inter¬ 
ferometers may be a better option. In the bottom panels of Fig. we show for such a 

possibility using 1, 10 and 100 pointings. Choosing ~ 80 m as the diameter for the array, the 
resolution is set at — 300. Here, the main design parameter is the field of view, FoV, which 
sets = 2n/s/FoV and is fixed by the effective size of each element, X / VF(^- For a 
‘dense array’, this is related to the number of elements, Ae ~ D\/d\^. Given that the maximum 
angular scale is set by the FoV, by adding more pointings, we simply diminish the signal variance 
by Ap, though the noise increases too, because the observation time Gbs —^ Gbs/A/- 

SKAl would use 15 m dishes, viz. £max — 60. Specifically, 254 dishes of SKAl-MID, 10,000 
hours and a survey area around 30,000 sq. deg. will correspond to ~ 2, as can be seen 

in Fig. H On the one hand, SKAl-MID will only go down to 350 MHz, whereas the minimum 
frequency assumed in Fig. ^ is 250 MHz. On the other hand. Camera et al. (2013) adopted a 
more conservative value for the system temperature, i.e. 30 K, whilst SKAl-MID will do better, 
with a system as cold as 20 K. Otherwise, SKAl-SUR will have a band from 900 to 350 MHz 
(0.58 < z < 3.06) and another from 1400 MHz to 700 MHz (0 < z < 1.1). Temperature will be 
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Figure 2: Forecast 1 a marginal error contours on ^ function of surveyed sky and total observation 

time, for a dish survey with Nd dishes (upper panels) and an interferometer making Np pointings (lower 
panels) (from Camera et al. 2013). 


slightly higher, at 50 K, for the former and 30 K for the latter. The effeetive number of dishes would 
be redshift dependent, most likely 64 x v for the firsthand with (z, v) = (0.6,30), (1.0,20), (2.0,9) 
and (3.0,5), whereas for the seeond band it would be 94 x v with (z, v) = (0.0,23), (0.5,10) and 
(1.0,6). For interferometrie surveys, we would need to wait for SKA phase 2, with the proposed 
‘aperture array’ system working below IGHz. It should be possible to achieve FoV ~ 1000deg^, 
thus reaching the cj(/^l) ^ 1 limit. See Santos et al. (2014) for a detailed description of SKA IM 
specifics. 


3.1.2 Primordial Non-Gaussianity Probed by the Bispectrum 

Compared to observable galaxies, HI is a weaklier biased tracer of the underlying matter dis¬ 
tribution, except for the largest scales where the scale-dependent correction due to PNG dominates. 
Therefore, on moderately large scales, we may model the HI bispectrum as the tree-level matter 
bispectrum modified by linear and non-linear bias factors, fini.t and bui ,2 ■ Here, we focus on the 
reduced HI bispectrum in redshift space, which is much less sensitive to other cosmological 
parameters (Sefusatti & Komatsu 2007). We have 


2 ?^ (^ 1 )^ 2 , ^3) 




AiP)] 




bm,2 

(^HI,l)^ 


(3.2) 


where is the reduced matter bispectrum predicted by the second-order perturbation theory, 
^o(jS) = 1 +2)3/3 -b j8^/5 and a® (j3) = 1 +2)3/3 + )3^/9, with jS = /(z)/fiHi,i the linear Kaiser 
factor. 

The relative importance of PNG in the HI bispectrum increases towards higher redshifts. This 
is very promising for measuring the primordial component from the HI bispectrum with SKAl- 
MID, which covers the whole redshift range from z ~ 0 to z ~ 3 (350 — 1420 MHz). In addition. 
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the large survey volume enabled by IM would result in a much better constraint on /nl- The full 
SKA, with its lower system temperature, better u-v coverage, and broader frequency band (hence 
broader redshift coverage and larger survey volume) will provide us with more stringent constraints 
on PNG than those achievable by other experiments such as the Tianlai project. 

For SKAl-MID, IM may be conducted with the dishes used individually in auto-correlation 
mode, being calibrated using interferometry. In this case, is limited by the Nyquist frequency, 
as well as the smallest scale above which we could trust the tree-level matter bispectrum. We 
adopt a non-linear scale cutoff, ^ni, for the tree-level matter bispectrum, by requiring the variance 
in the density contrast field at K/{2kn\) to equal 0.5 in each redshift bin. With Nd = 254 15 m 
diameter dishes, a survey area of 20,000 deg^, and a total integration time of 5,000 hr, we find 
^(/nl) ~ 45.7 and = 214.3 when we marginalise over fhe HI bias factors bun bui .2 

af each redshifl bin, while <7 (/nl) = 15.3 and <7 (/nl) = 61.8 if we assume a consfanf bias factors. 
Ofherwise, by using inlerferomelry wifh fhe full SKA, k^ax is sef by fhe non-linear scale k^u and 
we find much accurafe marginalised errors C7 (/^l) = 6.6 and <7 (/nl) ~ 55.4, or = 2.2 and 

o{Q) = 10.9. 


3.2 Intensity Mapping from the Epoch of Reionisation 

The PNG affects the clustering of the early star-forming galactic haloes responsible for cre¬ 
ating a network of ionised patches in the surrounding intergalactic medium during the EoR. This 
leaves a PNG imprint on the HI tomographic mapping in the intergalactic medium using its red- 
shifted 21 cm radiation. On large scales, where the typical size of ionised regions is much smaller 
than the scale of interest, we can neglect the non-linear effects of reionisation patchiness on the 21 
cm power spectrum. Then the 21 cm temperature power spectrum during the EoR can be written 
as 

P^^(k,z) = 5tIxI, [bui{k,z) + ^il]^P\k,z), (3.3) 

where dT b{z) = 23.88 / 0.02)y/0.15 -|-z)/10 mK, %i(z) is the global neutral hy¬ 

drogen fraction, and /Tk = k-n/k, i.e. the cosine of angle between the line-of-sight n and wave 
vector k of a given Eourier mode. 

The ionised density bias buw is the fundamental quantity derived from reionisation models, 
related to the neutral density bias bui by bui = (1 — Thii^hii) /Thi- The reionisation in the presence 
of PNG can be modelled using two independent methods as follows (D’Aloisio et al. 2013): 


1. Excursion-set model of reionisation (ESMR) - We can use a parameter ^esmr to characterize 
the efficiency of the local collapsed fraction of mass in luminous sources above some mass 
threshold in releasing ionising photons into the intergalactic medium (Eurlanetto et al. 2004). 
The full functions %i(z) and buii{k,z) are set by two parameters, i^esmr- 


2. Phenomenological model - Similar to the scale-dependent halo bias, D’Aloisio et al. (2013) 
derived a scale-dependent non-Gaussian correction to the ionised density bias, At?Hii(k,z), 
analogous to Eq. (2.2). It depends upon the scale-independent Gaussian ionised density bias 
^Hii(z)- Therefore, we can marginalise /nl over two phenomenological parameters, %i(z/) 
and bui\{zi), in each redshift bin Zi- 
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Both methods can be used to constrain PNG with the 21 cm power spectrum from the EoR. Mao 
et al. (2013) demonstrated that for a single frequency bin measurement, their constraints on 
are just as good. However, a model such as the ESMR can be used to combine multiple frequency 
bin measurements because it fixes the reionisation history (and therefore for a given value of 
/nl and (^esmr- Thus, it appears more promising. 

Eor a numerical evaluation, we assume that residual foregrounds can be neglected for > 
fc|| min = 27r/(yB), e.g. ^|| min = 0.055Mpc^^ at z = 10.1. The minimum k± is set by the minimum 
baseline, ^_L.min = 27rLmin/(ADa). We focus on the 21 cm signal on large scales ^max =0.15 Mpc^\ 
so that linearity conditions are met. We assume we can combine information from 7 frequency 
bins, each of 6 MHz bandwidth, covering the evolution of ionised fraction xuii ~ 0.06 — 0.70, 
corresponding to z ps 9.5 — 13.4 in the ESMR method. Using the Eisher matrix formalism, we 
forecast that the SKAl-EOW can constrain the local-type PNG, <7 (/nl) = 7.4, using the ESMR 
method. Eurthermore, if the sensitivity can be improved by 4 times better, then the full SKA can 
put as tight a constraint as <7 (/nl) = 4.7 by using the 21 cm power spectrum from the EoR alone. 


3.3 Multi-Tracer Technique 

Measuring clustering statistics on large scales has major the advantage of probing non-Gaussian 
effects on halo/galaxy bias where they are stronger. Unfortunately though, those same scales are 
limited by cosmic variance, i.e. the lack of enough independent measurements for the scales we are 
trying to probe, given the limited size of the volume that is observed. Seljak (2009) proposed a way 
to get around cosmic variance limitation by using different biased tracers of the underlying matter 
distribution. The main idea is that with at least two differently biased tracers, one can measure the 
ratio of these two biases to a level that is only limited by shot noise, hence beating cosmic variance. 
This is particularly sensitive when the bias of one tracer is much larger than that of the other. In 
order to understand this, let us assume that we measure two density maps at a given redshift, one 
for a biased tracer and another for the dark matter itself. Due to cosmic variance, there will be 
several cosmologies that are consistent with the dark matter map. However, the ratio of the two 
maps should give a direct measurement of the bias, with an uncertainty just given by the shot noise 
of the tracer (see also e.g. Abramo & Eeonard 2013). 

Eerramacho et al. (2014) developed this in the context of radio continuum surveys for SKA 
and its pathfinders. They explored the strong correlation that exists between radio galaxy types 
and their bias, which can be at first order linked to a single halo mass for each type. This seems 
to hold true for most types of radio galaxies, such as radio-quiet and radio-loud active galactic nu¬ 
clei, including ERl and ERll galaxies, and also star-forming galaxies, including starbursts. Using 
simulated catalogues, and taking into account all the statistical information from continuum galaxy 
surveys—namely the combined auto- and cross-correlation angular power spectra of multi-tracer 
galaxies—Eerramacho et al. (2014) showed that it is possible to strongly reduce the impact of cos¬ 
mic variance at large scales. By employing the Eisher matrix formalism, they forecast that this 
method will constrain the local-type PNG parameter up to an amazing accuracy of ~ 0.7, 

in the most optimal case where all the galaxy populations may be distinguished over the whole 
redshift range, and cj(/nl) ~ 2.9, in a much more conservative framework accounting for realistic 
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Figure 3: Constraints on obtained with the multi-tracer method as a function of the flux cut used to 
detect galaxies (from Ferramacho et al. 2014). The horizontal line represents the best constrain obtained by 
the Planck collaboration. 


limitations to the observational identification of each galaxy population.^ Note that the improve¬ 
ment obtained through the multi-tracer analysis is indeed significant, as the whole galaxy cata¬ 
logue without any galaxy type differentiation only allows for constraints on /nl with an error of 
^(/nl) ~ realistic scenario. Fig. presumes the forecast constraints as a function of flux 

cuf limif (see Ferramacho el al. 2014, for delails). 

3.4 Systematics Occurring on the Largest Scales 

The SKA should be able to provide a statistical detection of large scale effects due to its 
high sensitivity and large volumes probed, it will thus be crucial to control the systematics which 
will occur on these same scales. For the measurements we have been discussing, the following 
contaminants will need to be considered: 

• Masks - Several systematic effects related to partial sky coverage can result in biases and 
spurious large-scale signals. This is an issue for both galaxy surveys as well as IM. This 
effect can in principle be dealt with by means of inversion methods, as has been done by the 
CMB community for instance to address the masking of our own galaxy. 

• Foreground subtraction - For IM experiments, most foreground removal algorithms subtract 
out modes that are smoothly-varying in frequency, since this is how most galactic and extra- 
galactic foregrounds behave. In the process, some fraction of large scale cosmological modes 

^For more insight about the impact of redshift information for radio continuum surveys in the context of the SKA, 
see also Camera et al. (2012) and Raccanelli et al. (2014b), where in the latter work they used the ISW effect and redshift 
measurements from other experiments to constrain ~ 1. Comparable, encouraging results for SKA HI galaxy 

redshift surveys have recently been obtained by Camera et al. (2014c). 
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along the line of sight will also be subtracted. The effect should be most important for scales 
around the total bandwidth, but there will be a contamination on smaller scales too (Alonso 
et al. 2014). 

• Correlated noise - Auto-correlation observations are needed to access the largest angular 
scales with IM using the SKA. There is a noise term that has correlations in time, e.g. a non¬ 
white noise behaviour, with a drift that depends on the receiver design. This term cancels out 
in interferometric observations but is present in the auto-correlation signal from each dish. 
Therefore, strategies will have to be designed to remove this contribution. The situation can 
be particularly critical for large scale modes since this contamination shows up on time scales 
larger than the typical " 1/f" knee of the receiver, which can translate into striping artefacts 
on large angular scales. The usual solution is to design fast survey strategies capable of 
probing the required angles over a period that is shorter than the typical time-scale of the 
receiver. This will impose strong constraints on the receiver design that need to be accessed. 
Nevertheless other techniques should be explored. For instance the fact that the noise can 
also be a smooth function of frequency may allow for a solution in dealing with this effect. 

• Mosaicking - Large scale surveys will require several pointings (observations in different 
directions) in order to cover the required sky area—the so called ‘mosaicking’. This is true 
for any of the surveys, either continuum, HI galaxy or HI intensity mapping (the situation 
with SKAl-LOW is different due to its large field of view). Careful design of the scanning 
strategy needs to be made so that the noise is sufficiently uniform across the sky, in particular 
when overlapping the pointings. For example, fluctuations of this noise can generate differ¬ 
ent flux sensitivities across fhe sky and affecf fhe number of defecfed galaxies in a fhreshold 
survey. This could in furn generate spurious flucfuafions on angular scales (fhe sifuafion is 
less of a concern if we consider a high flux cuf in fhe galaxy survey, such as 10a). The same 
issue can arise if fhere are non-negligible flucfuafions on fhe calibrated signal. These prob¬ 
lems can in principle be deal! wifh fhrough fhe mosaicking approach as well as an accurafe 
calibrafion process. Moreover, fhese flucfuafions should be uncorrelafed on large angular 
scales which should allow fhem fo be distinguished from fhe signal. 

4. Conclusions 

The access fo observations on ulfra-large scales fhaf are possible wifh fhe SKA enables a range 
of unique science goals al fhe forefronl of cosmology. They include conslrainfs on DE and MG, 
amongsl fhe key science goals idenfified by fhe Commillee for a Decadal Survey of Asfronomy 
and Aslrophysics (National Research Council, 2010). The conslrainfs possible wifh fhe SKA will 
rival Ihose fhe Euclid salellile will deliver (Eaureijs el al. 2011; Amendola el al. 2013), bul will be 
complemenlary especially where fhe systematic effecls are concerned. The SKA will also be able 
fo use relalivislic effecls fo lesl GR on fhe largesl scales, and if will surpass fhe conslrainfs on PNG 
obfained by fhe Planck satellite. The SKA will probe fhe inflalionary dynamics lhal look place in 
fhe Universe jusl fractions of a second afler fhe Big Bang, and if will also lesl whelher gravity was 
modified during lhal epoch. 
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This Chapter, together with others in this review, elearly shows that the SKA will be able to 
probe the foundations of our eosmologieal models, the evolution of the Universe and the nature of 
its eonstituents. Moreover, as we have deseribed, this will be aeeessed by a range of SKA surveys 
whieh will probe the underlying signal in different ways. The SKA will eontribute in unique 
ways to broaden our understanding of the eosmos, and the SKA data will beeome a key tool for 
eosmology in the 21st eentury. 
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